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If the number of lattice sites is odd, a quantum particle hopping on a bipartite lattice with random 
hopping between the two sublattices only is guaranteed to have an eigenstate at zero energy. We 
show that the localization length of this eigenstate depends strongly on the boundaries of the lattice, 
and can take values anywhere between the mean free path and infinity. The same dependence on 
boundary conditions is seen in the conductance of such a lattice if it is connected to electron reservoirs 
via narrow leads. For any nonzero energy, the dependence on boundary conditions is removed for 
sufficiently large system sizes. 

PACS numbers: 72.15.Rn, 71.30.+h, 64.60.Fr, 05.40.-a 
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I. INTRODUCTION 



Zero-modes, wavefunctions at zero energy, often arise 
in problems when a quantum particle moves in _a 
background with a nontrivial topological structure E~lHI 
Quantum fluctuations associated to these zero-modes 
have dramatic physical consequences. They appear both, 
in field theories, related to chiral and parity anomalies ,u 
and in lattice regularization of field theoriesjj and have 
applications to a wide range of areas in physics: Chi- 
ral symmetry bpeaking in (l-l-l)-space-time quantum 
electrodynamics,0 edge states aloiiffi-|the boundary of a 
disk threaded by a magnetic flux,clEl singular contribu- 
tions to the Hall conductance from electrons hopping on 
a square lattice in the presence of a uniform magnetic 
fieldjLi superconductivity of a cosmic string,Ll localization 
of a fractkm al . phaijge at a domain wall in a charge den- 
sity wave 
rent in ^He-A 
semiconductors 
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induction of a persistent mass cur- 
antipfaase bojindaries in narrow-gap 
surfacell3 (edgeEj'ES) states in a super- 
condjuctor with dx2_y2 {-\-idxy) symmetry or in a chiral p- 
waveE2l supercatijductor, edge states in nanographite rib- 
bon junctions,L3 and itinerant-electron ferromagnetism 
in the repulsive Hubbard modelE3 

Related nontrivial topological structures can also exist 
in random matrix theorytil and in the problem of An- 
derson localization when the disorder possesses a spe- 



cial symmetry. In this context, an almost half a century 
old exarnple is that of a one-dimensional chain with link 
disorder.E3 Here, zero energy corresponds to the center 
of the energy band, and an eigenfunction with zero en- 
ergy is guaraeteed to exist if the number of sites in the 
chain is pdd.Ej In the thermodynamic .limit, the density 
of statefS and the localization lengthHa are singular 
at zero energy, wheTje^s .correlation functions of the lo- 



CTjpas .c 

Mm 



cal density of statesE3l£3 are algebraic functions. This 
anomalous behavior is an example of a strongly random 
critical point .lO Its origin is rooted in the stochastic prop- 
erties of the zero-modes supported by the Dirac equation 
in one-space dimension in the, backg round of a white- 
noise correlated random mass.E3c3l2a More recently, a 
two-dimensional random Dirac Hamiltonian with white- 
noise correlated U{\) x SU{N) randon^-weetor potential 
was shown to be critical at zero energy.E3c3 As with the 
stochastic model of Dyson, this critical behavior can be 
ascribed to the stochastic properties of zero-modes sup- 
ported by Dirac equations in two-space dimensions in 
the backgrjQJind of white noise correlated random vector 
potentials .El'OCJ 

While the critical behavior of the (continuum) Dirac 
equation is related to zero-modes in the infinite system, 
the existence of zero-modes for Dyson's stochastic model 
of a quantum particle hopping on a bipartite lattice with 
link disorder can also occur for a finite system size. The 
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Hamiltonian for this system is 



(1.1) 



where i and j label the lattice sites on a cartesian grid 
in d dimensions, say, and the hopping matrix elements 
tij are nonzero for nearest-neighbors only. Examples of 
bipartite lattices are depicted in Fig. |l]. In general, the 
tij will have a small random component in addition to 
an average t which sets the width of the spectrum of Ti,. 
In this paper we refer to Eq. (1.1) with the random tij 
as the "random hopping model" . The special case when 
it is only the phase of the hopping amplitude tij that is 
random is also known as the random flux problem. 

For the Hamiltonian (1.1), the existence of the zero- 
modes follows from the existence of a "sublattice" (or 
"chiral") symmetry. This symmetry follows when the 
lattice is divided into two sublattices A and B such that 
the hopping matrix tij only connects sites from the two 
sublattices, but not sites of the same sublattice. For the 
example of Fig. ||, the sublattices A and B correspond to 
the white and black sites, respectively. In a matrix form 
and after a relabeling of indices, the eigenvalue problem 
TYIV") = elV') can be rewritten as 



t 



"AB 



AB 




^A 
i'B 



(1.2) 



where V'A and iJjb denote the wavefunction on the lattice 
sites of the sublattices A and i?, respectively. Then, de- 
noting the number of sites in the sublattice s A and B by 
Na and Nb, counting dimensions in Eq. (1.2) immedi- 
ately yields that thejaiimber of linearly independent zero- 
modes IS \N A- Nb\ -Ell To see this, note that if Na> Nb, 
ipA obeys an underdetermined set of linear equations, 
while ipB obeys an overdetermined set of equations. (For 
all lattices shown in Fig. ^ there is one zero-mode with 
support on sublattice A.) 

Addition or removal of a single site changes the num- 
ber of zero-modes, as it changes the difference \Na — Nb\ 
between the numbers of sites in the two sublattices. At 
the same time, the singular behavior of the density of 
states and the localization length in the random hopping 
model in the thermodynamic limit are considered "in- 
trinsic" properties, i.e., they are derived from continuum 
models and should not depend on boundary conditions 
of the lattice. Hence, whilst both the existence of zero- 
modes for lattices with boundaries and the singular be- 
havior of the localization length are manifestations of the 
same sublattice symmetry, they are so in very different 
ways. One might even ask to what extent the localization 
length of zero-modes is representative for the "intrinsic" 
localization length of the random hopping problem or 
random flux problem on a lattice without boundaries. 
This is the question addressed in this paper. 

Our answer is that the localization length of the zero- 
modes for lattices with boundaries is not an "intrinsic" 
property of the random hopping model. After a brief re- 
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FIG. 1: Three examples of a lattice with different boundary 
conditions. The sublattices A and B correspond to the white 
and black sites, respectively. In all three examples, the num- 
ber Na of sites on sublattice A is one more than the number 
Nb of sites on sublattice B. (a) Conventional rectangular 
shaped wire with N — 3 and odd L ^ N. (b) Boundary 
conditions with respect to (a) have been changed by adding a 
white site to the left and adding a black site to the right, (c) 
Boundary conditions with respect to (a) have been changed 
by adding a pair of white and black sites to the left. 



support this conclusion in Subsec. [IB by analytical and 



numerical solution of the problem in the case of a wire 
geometry: Depending on the boundary conditions, zero- 
modes exhibit a range of localization lengths, the smallest 
one being of the order of the mean free path. In Subsec. 
11 C we then show that this extreme sensitivity to bound- 
ary conditions is an anomaly corresponding to the special 
case £ = 0. For any energy e ^ there is a unique local- 
ization length if the system size is sufficiently large. For 
sufficiently small e, this unique localization length coin- 
cides with the largest of the possible localization lengths 
at e = 0. We discuss higher dimensional examples in Sec. 
HI and conclude in Sec. 



view of the transfer matrix formalism in Subsec. II A, we 



II. QUASI-ONE-DIMENSIONAL GEOMETRY 

In this section we consider a two-dimensional lattice, 
N sites wide and L ^ N sites long. Examples for = 3 
are given in Fig. |^ whereby sites on sublattices A and B 
are colored in white and black, respectively. The number 
of zero- modes equals \Na — Nb\- For definiteness, we 
assume A^^i > Nb, so that all zero- modes have support 
on sublattice A. 
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quire a dependence on the position y along the wire. The 
relation between the a„(?/) at different positions along the 
wire can be expressed through the transfer matrix M. , 
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'(y) = X! Mnm{y,y')am{y')- 



(2.1) 
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FIG. 2; Disordered quantum wire of even length L (dotted 
region) coupled to left and right reservoirs through leads of 
the same width as the wire. We have chosen leads without 
transverse hopping for technical convenience. 



In this section, we will compare zero-modes in wire 
geometries with different boundary conditions at the two 
ends of the wire, as is shown, e.g., in Fig. |l|. We will 
establish that the localization length of these zero-modes 
cannot be thought of as being intrinsic, i.e., independent 
of the boundary conditions even as the thermodynamic 
limit L — > oo is taken. We will then turn our attention to 
finite energies and, supported by a numerical solution of 
the problem, argue that an intrinsic localization length 
at arbitrarily small but finite energies does indeed exist. 
The order in which the limits L — > cxd and e ^ are 
taken is thus essential for the extraction of an intrinsic 
localization length at zero energy. 

The sublattice symmetry singles out the band center 
in that, under an appropriate choice of boundary condi- 
tions, an exact energy eigenfunction at that energy can 
be constructed for any realization of the disorder. This 
is not true of any finite energy e in a closed and finite 
system. Therefore we will proceed in two steps. We first 
compare the localization length of zero-modes in a closed 
system with the exponential decay lengths for transmis- 
sion probability of plane waves in an open system and es- 
tablish that they are equal. Then we study how the trans- 
mission probability of plane waves is changed when the 
energy becomes finite. In both steps we use the transfer 
matrix formalism, which can deal with open and closed 
systems in a unified way. 



A. Transfer matrix formalism in an unbounded 



1. Plane-wave representation 

In the absence of disor der, the eigenfunctions of the 
Schrodinger equation (l.f ) in a wire geometry as depicted 
in Fig. H are plane waves. At zero energy, there are TV 
independent wavefunctions for plane waves traveling to 
the right, and N independent wavefunctions traveling to 
the left. An arbitrary wavefunction can be expanded in 
the basis of plane waves. In the presence of disorder, the 
plane- wave expansion coefficients a„, n = f, . . . ,2N ac- 



Current conservation and the siiblattice symmetry imply 
that Ai can be parametrized asC3 

U \ / coshX sinhX \ / \ 
M=\ . (2.2) 

U J \ sinhX coshX J \ V J 

The 2x2 grading displayed here is that of right and 
left moving plane waves. The N x N matrices U and V 
are unitary (/? = 2), symplectic (/? = 4), and orthogonal 
{(3 — 1) when the hopping amplitudes are complex, real 
quaternions, and real valued, respectively. The N x N 
matrix X = diag(xi , . . . , xn) is real valued and diagonal. 
In the absence of the sublattice sympietry, the transfer 
matrix has a similar parametrization.c3 The main differ- 
ence between the cases with and without sublattice sym- 
metry is that without sublattice symmetry the Xji can 
always be chosen positive, while with sublattice symme- 
try both positive and negative Xn appear. 

The localization properties of the disordered wire are 
encoded in the transfer matrix M. They are dominated 
by the diagonal matrix X in Eq. (2.2). The distribution 
of the Xn at zero energy has been studied in Ref. ^ for 
the case of a wire with random hopping only. In terms 
of the mean free path I of the wire, it was found that 



L 



n — 1, . 



(2.3) 



if L ^ Nl, up to fluctuations of relative order 
{Nl/ L)^^/"^ ^ where the Lyapunov exponents f/^„ (the 
inverse localization lengths) are given by (to leading and 
subleading order in iV, see Ref. Eq) 



f [3{N + 1 - 2n) 



71 = 1, . 



,N. 



^„ (/3iV + 2 - p)l ' 
For an infinite wire, the conductance is given by 



(2.4) 



2e2 

G = —— cosh " X. 



2 „ — 



h ^ 

n=l 



2e^ 

^9- (2.5) 



Hence, for an infinite wire it is only the smallest in mag- 
nitude of the Lyapunov exponents that governs the ex- 
ponential decay of the conductance. Its inverse is thus 
identified with the localization length ^ of the system. 



iPN + 2 - P)l/f3 for N even. 



for N odd. 



(2.6) 



As we will find below that all localization lengths ^„ can 
serve as localization lengths for zero-modes in a finite 
sized wire, we will refer to their maximum ^ as the local- 
ization length for an infinite wire. 
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FIG. 3: Choice of the vectors of Eq. (2.7). 



2. Site representation 

An alternative representation for the transfer matrix is 
obtained using a site representation for the wavefunction 
instead of an expansion in plane waves. The sublattice 
symmetry becomes manifest in the site representation if 
the wavefunction elements are arranged in A'^-component 
vectors ipy containing elements of one sublattice only. In 
that notation, the Schrodinger equation (1.1) reads 



(2.7) 



The index y labels the coordinate along the wire and the 
iV-component vectors ipy contain wavefunction elements 
on sites of sublattice A if j/ is odd and of sublattice B if 
y is even. A possible choice for the vectors V-'y is shown 
in Fig. ^. The N x N matrices Ty contain the hopping 
amplitudes between adj ace nt sites. 

The solution of Eq. ( ^.7|) can be represented in terms 
of a transfer matrix as well, 



Miy,y') 



M{y,y') 



n 



i^y, 



1 



-'-m-l 

1 



(2.8) 



We note that, if y — > is even, the transfer matrix 
Ai{y,y') is block diagonal at zero energy. In that case, 
Mly, y') can be parametrized as 



Miy,y') = 



Ue^V 
T-^lUe 



(2.9) 



Here U and V are orthogonal (unitary, symplectic) ma- 
trices ioT (3 — 1 {(3 = 2, 4) a nd X is a diagonal matrix. 
The grading used in Eq. ( |2.9| ) corresponds to the division 
into sublattices A and B. As the transfer matrices of Eqs. 



( p.8[ ) and (2.1) are related by a simple basis transforma- 
tion, the distribution of the eigenvalues Xn of the matrix 
X is also given by Eqs. (2.3) and ( ^.4[ ) above. Hence, 
as long as we are interested in the Lyapunov exponents 
only , we can choose freely between the site representation 
(p.8[) and the plane wave representation (|2.2|). 



B. Bounded wires: zero energy 

1. Wavef unctions 

Zero-modes are solutions to Schrodinger equation ( |l.2| ) 
at zero energy. To see how they are constructed, let us 
first inspect the case of Fig. |l|(b) in detail. In this case, 
the zero mode is a wavefunction with support on the 
white sites only. To construct it, we start with three ini- 
tial values (seeds) for wavefunction on the three leftmost 
white sites. We now construct the wavefunction on all 
other sites in the same wa y as one finds the transfer ma- 
trix in site representation (2.8): By applying Schrodinger 
equation to the leftmost black site of the middle row one 
obtains the value of the wavefunction on the penultimate 
leftmost white site of the middle row. Repeating this 
process column by column for all black sites that have a 
white site to its right, we can construct a wavefunction 
supported on all white sites. For it to be a zero-mode, 
the Schrodinger equation must also be satisfied on the 
two rightmost black sites that were not used to propa- 
gate the wavefunction on the white sites. For Fig. wOo) 
these are the two rightmost black sites. Application of 
Schrodinger equation on each of those sites yields two 
linear constraints for the wavefunction elements on the 
rightmost white sites. Both constraints can be satisfied 
since they are implemented linearly on the three seeds of 
the wavefunction on the left end. 

This example suggests a three-step recursive method 
to obtain zero-modes. First, the zero-mode has support 
on one sublattice only. Second, N independent numbers 
that make up the A^-component vector ipAL are assigned 
to the values taken by the wavefunction on the leftmost 
sites of each row that belong to sublattice A. Solution of 
the Schrodinger equation on all sites of sublattice B ex- 
cept for those without a white nearest-neighbor to their 
left or right (i.e., except for sites of sublattice B that are 
at the left or right ends of a row) allows to propagate 
recursively the wavefunction to the right column by col- 
umn. The rightmost N values of the wavefunction thus 
constructed build the vector ipAR- The relation between 
the vectors ipAL and ipAR that is thus obtained can be 
expressed as 



AB. 



'AL- 



(2.10) 



The N X N matrix M is nothing but the counterpart to 
the upper left block of the transfer matrix in site repre- 
sentation. Hence, by Eq. (2.9), it has the polar decom- 
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position 



(2.11) 



where U and V are orthogonal (unitary, symplectic) ma- 
trices for (3 = 1 {P — 2,4), and X is a diagonal matrix 
with eigen valu es x„ whose statistics are given by Eqs. 
(U) and (U) above. 

Third, in order to have a true zero-mode, the 
Schrodinger equation must be obeyed on the remaining 
Nb l and Nb r black sites that do not have white nearest- 
neighbors to their left and right, respectively. This gives 
Nbl independent constraints to be satisfied by the el- 
ements of il^AL and Nbr independent cons traint s to be 
satisfied by the elements of i^ar- As Eq. ( 2.10 ) allows 
for N independent solutions, the number of independent 
zero-modes is equal to 

# zero- modes ^ Na - Nb = N - {Nbl + Nbr) ■ 

(2.12) 

The equality Na- Nb = N - (Nbl + Nbr) followed 
since the lattice topology is such that the only black sites 
without left or right white nearest-neighbors are located 
at the far left and far right ends o f the wire, respectively. 
We emphasize that the criterion ( ^.12| ) for the existence 
of zero-modes does not rely on the quasi-one-dimensional 
assumption N <^ L or on the assumption that there be 
only one "tran svers e" direction. The only ingredients 
needed for Eq. ( 2.12 ) to hold are the boundary conditions 
at the end of the disordered region, i.e., the topology of 
the "wire" . The simplifying feature brought by the quasi- 
one-dimensional limit N/L <C 1, is that the statistical 
properties of the transfer matrix M in Eq. ( 2.10| ) are 
known. We will exploit this knowledge below. 

How do the wavefunction elements at the left and right 
ends of the wire compare? To answer this question, we 
first look at the geometric mean </> of ipAL and "0^1^?: 



AR- 



(2.13) 

In terms of this geometric mean, the Nbl constraints on 
ipAL and the Nbr constraints on ipAR. can be written in 
the form 



Cbl e 
Cbr e 



'X/2 
-X/2 



= 0, 



(2.14) 



where Cbl and Cbr are {Nbl x N) and {Nbr x N) 
matrices with coefficients of order unity, respectively. In 
the localized regime L ^ Nl, the spaced by an 

amount of order L/Nl ^ 1, so that the coefficients in 



simp 



Eq. (2.14) differ considerably in magnitude. To see what 



ifications this brings about, we look at the first row 



of Eq. ( |2T4| ), 



(CBL)iie-"^/'(/.i + . . . + {CBL)lNe-^"/^^N = 0. 



(2.15) 



According to Eqs. (2.3, 2.4), the coefficient 
{CBL)ine~'^"^^ is a random mmiber that ffuctuates 
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FIG. 4: Example of a disordered quantum wire of even length 
L (dotted region) coupled to left and right reservoirs via point 
contacts. In this example the number of left end points on 
sublattice A is Nal = 2 and the number of right end points 
on sublattice B is Nbr ~ 2. The number of right end points 
Nar on sublattice A vanishes as does the number of left end 
points Nbl on sublattice B. If L is chosen to be odd, Nal = 
Nar = 2 whereas Nbl ~ Nbr = 0. When L is odd, a 
zero-mode can only be supported on sublattice A and the 
conductance must vanish identically. 



around e^'^/^^^"-'. Since by Eq. ( p^ ) xn is smaller than 
all other a;„ by an amount of at least L/N l » 1, we thus 
find that the left hand side of Eq. (2.15) is dominated 
by the last term, so that we conclude 0jv — 0, with 
exponential accuracy. Extending this argument to the 
first Nbl and the last Nbr rows of Eq. ( 2.14 ) we infer 
that 



0, 
0, 



n 
n 



N -Nbl + 1, 
1, • • ■ , Nbr, 



,N, 



(2.16) 



( p. 13 ), the only localization lengths available tollie 
modes are ^„ with n ~ Nbr -t- 1, . . . ,N — Nbl^ ] 



respectively, again to exponential accuracy. Conversely, 
the only nonzero elements of (f) are (j)n with n = Nbr + 
1, . . . ,N — Nbl, to exponential accuracy, so that, by Eq. 

1e zero- 
In our 

notations, a negative localization length corresponds to a 
wavefunction exponentially localized near the left end of 
the wire (since then \^pAL\ ^ iV'Afil in that case), while 
a positive localization lengths correspond to a wavefunc- 
tion exponentially localized near the right end of the wire. 
A divergent localization length (which can occur for odd 
N) signals a zero-mode that is critical (|V'al| and \iPar\ 
comparable in magnitude to exponential accuracy). 

We have verified this scenario by numerical implemen- 
tation of above recursive construction of zero-modes in 
geometries depicted in Fig. flifor various choices of N 
and of boundary conditions.EsI The agreement found is 
excellent. 



2. Transmission probability 

A different method to probe the effect of boundary con- 
ditions on the localization length of the random hopping 
model is via the transmission probability of a lattice cou- 
pled to electron reservoirs via ideal leads. Boundary con- 
ditions play a role once the width of the leads is smaller 
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FIG. 5: Four difFerent bounded wires N — 5 sites wide 
and L sites long. In all cases L is chosen even and the wire 
is connected to the reservoirs by single-channel leads. The 
boundary conditions are specified by: (a) Nal = Nar = 
Nbl = Nbr = 2. (b) Nal = Nbr = 1, Nbl = Nar = 3. 
(c) Nal = Nbr = 0, Nbl = Nar = 4. (d) Nal = 0, 
Nar = 3, Nbl = 4, Nbr = 1. 



than the width of the sample lattice, as is shown, e.g., in 
Figs. |, I and |. 

With ideal leads attached to the left and right ends, 
zero-modes with support on sublattice A can coexist with 
zero-modes with support on sublattice B. In fact, since 
a traveling wave at zero energy has support on both sub- 
lattices, conductance through the sample is only possible 
if both types of zero- modes exist Jla Using the same ar- 
guments as for the zero-modes in a closed system, the 



FIG. 6: Three difFerent bounded wires = 6 sites wide 
and L sites long. In all cases L is chosen even and the wire 
is connected to the reservoirs by single-channel leads. The 
boundary conditions are specified by: (a) Nal ~ Nbr = 2, 
Nbl = Nar = 3. (b) Nal = Nbr = 1, Nbl = Nar = 4. 
(c) Nal = Nbr = 0, Nbl = Nar = 5. 



possible localization lengths for zero-modes in the pres- 
ence of leads can be found from counting the number of 
end points belonging to each sublattice. More precisely, 
let Nb l and Nbr be the number of sites of sublattice B 
without a white nearest-neighbor to their left and right, 
respectively. For all geometries under consideration, Nb l 
(Nbr) match the number of end points on the left (right) 
end of the wire that belong to sublattice B (sites con- 
nected to leads are excluded here). Now, the quasi- 
one-dimensional localization lengths of the zero-modes 
with support on sublattice A are Cnbr+i^ ■ ■ ■ i^n-Nbl- 
Similarly, the available quasi-one-dimensional localiza- 
tion lengths for zero-modes with support on sublattice 
B are Cnar+i^ ■ ■ • ,^n-Nal^ where Nal (Nar) denotes 
the numbers of end points on the left (right) end of the 
wire that belong to sublattice A (again, sites connected 
to leads are here excluded). Whether or not the quasi- 
one-dimensional limit apphes, if Nbl + Nbr > N ot 
Nal + Nar > N, there are no zero-modes with support 
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TABLE I: Maximum localization lengths and ^g, of zero- 
modes on sublattices A and B, respectively, for the four dif- 
ferent geometries depicted in Fig. H. The minimum of and 
denoted controls the conductance. The entry = 
for geometry (d) indicates that no zero-mode is supported on 
sublattice A. Correspondingly, the conductance vanishes at 
zero energy in this geometry as is implied by the vanishing 
localization length f * — 0. 
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(d) 
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TABLE IL Maximum localization lengths ^\ and ^g, of zero- 
modes on sublattices A and B, respectively, for the three dif- 
ferent geometries depicted in Fig. H. The minimum of and 
^B, denoted controls the conductance. 
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on sublattice A or B, respectively, and hence no traveling 
waves and no conductance through the sample. If both 
Nbl+ Nbr < N and Nal + Nar < N, there is a finite 
conductance through the sample. In that case, the quasi- 
one-dimensional exponential decay length of the conduc- 
tance, denoted ^* , is the shorter one of the maximum of 
the decay lengths for zero-modes on the two sublattices. 
We give in Tables || andj| the values taken by ^* for the 
four geometries of Fig. ^and the three geometries of Fig. 
^, respectively. 

Thus, we find that the same range of localization 
lengths shows up in the exponential decay of wavefunc- 
tions in a closed system and of the conductance G in an 
open system, when the wire is coupled to the electron 
reservoirs via point contacts. This is in stark contrast 
to the case of an "infinite" quantum wire (i.e., a wire 
without point contacts at both ends as depicted in Fig. 
H), where only the largest localization length determines 
the conductance. It is also in contrast to the case of 
a quantum wire with standard diagonal disorder, where 
the boundary conditions have no effect on the exponen- 
tial decay length of the conductance. 

Again, we have verified this scenario and found excel- 
lent agreement between numerics and theoretical expec- 
tations. 



C. Bounded wires: nonzero energy 

To study the importance of boundary conditions at 
finite energy, we have calculated numerically the conduc- 
tance of a lattice with random hopping amplitudes and 



point contacts as a function, nf energy using the recur- 
sive Green function method.EBta In our numerical sim- 
ulations we chose real-valued nearest-neighbor hopping 
amplitudes tij in the disordered region uniformly and in- 
dependently in the intervals 1—S<tij < l+S for hopping 
in the longitudinal direction and t{l ~ S) < tij < t{l + S) 
in the transverse direction, with S = 0.2 and t = 0.6. 
With this choice the mean free path I is about 65 lattice 
spacings, as can be estimated from a fit to (In 5) versus 
L for large L and large energy e (« 10~^). In the perfect 
leads we take (5 = and t = 0. The numerical data are 
obtained after averaging over 10^ samples. The size of 
error bars for {\ng) is estimated to be less than 1%. A 
more detailed account of our implementation of the re- 
cursive Green function method can be found in Refs. |5^ 
and ^ The disorder strength in the numerical calcula- 
tions presented here is chosen the same as in Ref. so 
that a comparison of the results is possible. 

We have calculated (In g) versus L for four different 
boundary conditions for a quantum wire of width N = 5 
and three different boundary conditions for a quantum 
wire of width = 6. The boundary conditions are shown 
in Figs. H and ^ Our results are shown in Figs. |7[ ^, ||, 
and 1^. 

Upon increasing the energy away from the band center 
£ = 0, a crossover from the chiral (orthogonal) symme- 
try class to the standard (orthogonal) symmetry class is 
expected to take place. For wires without point contacts 
at the left and right ends, this crossover was studied by 
three of the authors in Ref. There, we found that the 
crossover to the standard orthogonal class happens for 
e ~ Ec, where 



(2.17) 



is the Thouless energy for a localization volume of lin- 
ear dimension ~ Nl. (The relevant localization length is 
the smaller one of the localization lengths in the chiral 
and standard symmetry classes) . For our calculations we 
estimate that, with vp ~ 2, Ec ~ 10"'^ for N = 5 and 
A^ = 6. The largest energy we consider in our calcula- 
tions, e = lO^^is well inside the standard orthogonal 
class, see Fig. |l^. For that largest energy, the four or 
three curves of (In g) versus L/l in Figs. ^ and ||, respec- 
tively, that correspond to the different boundary condi- 
tions are indistinguishable. The same conclusion can be 
reached from Figs. ^ from which one can infer the Lya- 
punov exponents 



1 d{\ng) 

2 dL 



L:»Nl, 



(2.18) 



or from Figs. |T^, where we showed the energy dependence 
of (In 5) at a fixed length. This agrees with the con- 
ventional understanding that, in the absence of the sub- 
lattice symmetry, the localization length is an intrinsic 
property of the wire, and hence boundary-independent. 

The two other energies we considered (e = 10~^ and 
10~^) are both much smaller than Sc, i.e., well inside the 
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FIG. 7: Dependence on the wire length L, L always even, of (Ing) for a quantum wire with random hopping only, connected 
to the reservoirs via single-channel ideal leads, and with the boundary conditions specified in Fig. |^. The quantum wire is 
N = 5 sites wide. The curves represent different values of the energy e and different choices for the boundary conditions at the 
left and right ends of the wire, as depicted in Fig. pi 



chiral class. For short lengths, we see the same strong 
dependence on boundary conditions that was predicted 
for the zero-modes in the previous subsections. A quan- 
titative verification of the predictions is found from Fig. 
^, where the rescaled Lyapunov exponents 

{N + 1)1 d{h-ig) 
2 dL~ 

take even integer values for odd N and odd integer values 
for even N fo r short L (but still L ^ Nl), in agreement 
with Eq. (2.4) with (3 — 1. For larger lengths, however, 
the dependence on the boundary conditions is lifted, and 
the Lyapunov exponents are the same for all boundary 
conditions considered as is illustrated in Fig. ^ at energy 
e = 10~^. For a sufficiently large length of the wire and 



upon decreasing the energy, —{N + l)ld{lng)/2dL ap- 
proaches for 5 and 1 for A = 6, irrespective of 
the boundary conditions. Again, this is well illustrated 
by Fig. ^ at energy e — 10~^. Alternatively, for N — 5 
(odd A^) and energy e = 10~^, (Ing) in Figs. 0(b-c) is 
found to depend linearly on length for small L with a 
boundary-condition dependent slope, while (In 17) has a 
curvature consistent with a (L/Niy/"^ dependence, as is 
appropriate for critical conductance statistics. For en- 
ergy e — 10~^ the L-dependence of (Ing) is linear for 
large L, but with a localization length that is signif- 
icantly larger than for e = 10^^. Such an enhanced 
localization length is characteristic of the crossover be- 
tween the chiral andrStandard classes for a quantum wire 
without boundaries.^ For A = 6 (even A^) and energies 
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FIG. 8: Dependence on the wire length L, L always even, 
of (In g) for a quantum wire with random hopping only, con- 
nected to the reservoirs via single-channel ideal leads, and 
with the boundary conditions specified in Fig. ^ The quan- 
tum wire is A'^ = 6 sites wide. The curves represent different 
values of the energy e and different choices for the boundary 
conditions at the left and right ends of the wire as depicted 
in Fig. |. 



e = 10~^, 10~^, (Ing) in Figs, [^(a-c) decreases linearly 
with length, but with different slopes for small and large 
L. These slopes correspond to exponential localization 
controlled by a boundary-condition dependent zero mode 
and to exponential localization in the chiral orthogonal 
class (or, strictly speaking, the crossover between the chi- 
ral and standard orthogonal classes) in an infinite wire, 
respectively. (Note, however, the large energy-dependent 
crossover lengths and the non-monotonous length depen- 
dence of the Lyapunov exponents at intermediate length 
scales in Figs. 0-lO[) Hence, from the numerical cal- 
culations we conclude that for a finite energy, the Lya- 
punov exponents lose their dependence on the boundary 
conditions if the wire is sufficiently long. The typical 
conductance itself, exp((ln(7)), retains a strong depen- 
dence on the boundary conditions for sufficiently long 
wires through its exponential prefactor as is illustrated 
by Figs. 1^ and However, the slope of (In 5) as a func- 
tion of L/Nl in the regime L/Nl ^ 1 is independent of 
the boundary conditions. 

The length scales where the Lyapunov exponents be- 
gin to cross over from the boundary-condition-dependent 
value characteristic of the zero-modes to the "intrin- 
sic" (smallest) Lyapunov exponent can be estimated as 
the length scales where the Thouless energy eTh{L) = 
g{L)A{L) is equal to the energy e, with g{L) ~ 
exp(— 2L/|^„|) the typical boundary-dependent dimen- 
sionless conductance of the wire {n — 1,... ,N) and 
A(L) = the mean level spacing of a wire with length 
L. Hence, 



l^nlln 



1,...,N. (2.19) 



No useful crossover length can be defined for the bound- 
ary condition of Fig. ra(d), where the zero-energy con- 
ductance of the wire is zero [see Table I and the down- 
turn of the traces with energies e = 10 , 10^^ for suffi- 
ciently small wire lengths in Fig. 1(d)]. Equation (|2T9| ) 
implies that the shorter |i^„| disappear at shorter wire 
lengths than the larger |^„|. All dependence on bound- 
ary conditions is removed and only the smallest of the 
Lyapunov exponents survives for lengths larger than the 



i.e 

E L 



beyond 



s,{N- 



-1)12 for N even. This 



second largest of the L^,, 
for N odd and beyond 

is well illustrated in Figs. ^, where the shortest local- 
ization length (corresponding to the highest Lyapunov 
exponent) disappears first, to be followed by the second- 
shortest localization length at a wire size that is roughly 
a factor 2 (N = 5) or 5/3 (N = 6) larger (see Tables || 
and respectively) . 

To summarize: While the exponential decay length of 
the conductance depends on the boundary conditions for 
e = even in the limit L — > 00, the exponential decay of 
the conductance is governed by the "intrinsic" (largest) 
localization length ^ for any finite energy e different from 
zero. In this sense, the remarkable dependence of the 
zero-mode localization lengths and of the zero energy 
conductance on boundary conditions can be considered 
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FIG. 9: Lyapunov exponents d(ln g) / dL for the curves shown 
in Figs. 1^ (upper panel) and ^ (lower panel). 



as an anomaly belonging to the case of energy being ex- 
actly equal to zero, not as something representative of 
the thermodynamic limit of the random hopping model. 



III. HIGHER DIMENSIONAL EXAMPLES 

The examples we have discussed so far pertain to a 
quasi-one-dimensional geometry. We would like to close 
with examples of lattices that are extended in two or 
more dimensions. For the three lattices shown in Fig. |l^, 
three different scenarios apply to the zero-energy conduc- 
tance between the left and right leads. 

For the lattice of Fig. |l^(a), the conductance decays 
exponentially with a decay length of the order of the 
mean free path of the system as can be seen by direct 
simulation of the conductance or by constructing recur- 
sively zero-modes on each sublattices. In the latter case, 



FIG. 10: Crossover as a function of energy for a fixed length 
L — 201 of the wire of {Ing) for N = 5 (upper panel) and 
= 6 (lower panel). All traces saturate to different finite 
values at very low energies except for trace (d) of the upper 
panel for which the conductance vanishes at the band center. 



a seed value is assigned first to the leftmost (rightmost) 
black (white) site having a pair of white (black) nearest- 
neighbor sites in the direction orthogonal to the leads. 
The Schrodinger equation is then solved on the white 
(black) sites moving to the right (left) column by col- 
umn. The zero-mode on black sites is exponentially lo- 
calized near the right contact and the zero-mode on white 
sites is exponentially localized near the left contact with 
a localization length of the order of the mean free path. 
The exponentially small conductance at the band cen- 
ter is found in all "diamond" -like structures of the form 
of Fig. pT|(a), irrespective of dimensionality and disorder 
strength. 

One possible generalization of the "diamond" -like ge- 
ometry of Fig. ^(a) is to couple sites on the lower left 
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FIG. 11: Three examples of two-dimensional lattices that 
have different size-dependencies of the conductance, see text. 



and upper right edges of the diamond to the reservoirs. In 
this case, the number of zero-modes scales with the width 
of the leads. To each seed on one edge corresponds a zero- 
mode exponentially localized about its vis a vis on the 
opposite edge. Linear superposition of these zero-modes 
yield traveling waves on the edges that are exponentially 
localized in the direction orthogonal to the edges. This is 
reminiscent of the Callan-Harvey effectE in field theory. 

For the "square" lattice of Fig. ^(c), the conductance 
is zero, as this geometry does not allow a zero-mo de on 
the black sites according to the arguments of Subsec. 11 B . 
(Recall that the condition of quasi-one-dimensionality is 
not needed to establish the conditions for the existence 
of zero-modes.) 

Finally, the lattice of Fig. ^(b) has zero-modes on 
white sites and on black sites that are believed to be 
critical (i.e., not exponentially localized). This is not 
surprising in a quasi-one-dimensional geometry for which 



the lattice only extends in the longitudi nal d irection, as 
can be seen using the arguments of Sec. IIB . More sur- 
prising is that the critical nature of the zero-mode seems 
not to depend on the transversal extension. Verges in 
Ref. ^ has studied numericall y t he conductance distri- 
bution in the geometry of Fig. |ll|(b) for a square lattice 
made up of up to 799 x 799 sites. His conclusion is that 
the probability distribution in the geometry of Fig. |ll|(b) 
is best fitted by the conductance distribution of a thick 
quantum wire with an odd number N of channels derived 
in Ref. ^ He thus concludes that the critical zero-mode 
in the two-dimensional geometry of Fig. [TT|(b) is quasi- 
one-dimensional in nature. 



IV. CONCLUSIONS 

In this paper we have investigated the dependence on 
boundary conditions of localization properties of the ran- 
dom hopping problem at the band center and its vicinity. 
At finite energies, localization properties are intrinsic, 
i.e., independent of boundary conditions in the thermo- 
dynamic limit. Remarkably, this is not true anymore 
precisely at the band center where both the transmission 
probability of a disordered region connected to reservoirs 
by single-channel leads and the spatial decay of zero- 
modes in closed systems are highly sensitive to the choice 
of boundary conditions even as the thermodynamic limit 
is taken. This sensitivity to boundary conditions was 
quantified analytically in quasi-one-dimensional geome- 
tries. In particular, the conditions under which zero- 
modes are critical were given. In higher dimensions, one 
must rely on numerical simulations to study localization 
properties of zero-modes. However, the task is some- 
how simplified by an explicit recursive construction of 
zero-modes that we showed applies to a large class of ge- 
ometries. It is an interesting open problem to determine 
conditions for criticality of zero-modes beyond quasi-one- 
dimcnsional geometries and whether a field theoretical 
description of intrinsic critical properties in terms of crit- 
ical zero-modes in the spirit of Refs. and ^ ap- 
plies to the two-dimensional random hopping problem. 
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